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By Stephen H. Maslen 


SUMMARY 

The first-order solution for the laminar compressible boundary- layer 
flow over a flat plate at constant wall temperature is given. The 
effect of slip at the wall as well as the interaction between the 
boundary- layer flow and the outer stream flow are taken into considera- 
tion. The solution is obtained explicitly in terms of the known zero 
order, or continuum, solution. No assumptions regarding the Prandtl num- 
ber or viscosity- temperature lav need be made. It is found that the 
first-order solution gives a decrease in heat transfer and, for super- 
sonic flow, an increase in skin friction. For subsonic flow there is no 
first-order shear effect. The change in heat transfer is due to slip 
and the change in friction is due to the interaction of the zero- and 
first-order velocities at the outer edge of the boundary layer. 


INTRODUCTION 

With very high-altitude high-speed motion becoming of practical 
interest, the behavior of air flow in rarefied and semirarefied gases 
becomes of great importance. In this connection, as has been discussed 
by Tsien (reference l), one can define four regimes of fluid flow. 

These may be termed, in order of increasing mean free path of the fluid 
molecule: continuum, slip, intermediate, and free molecule flows. Con- 

tinuum flow, where the mean free path of the fluid is negligible com- 
pared with the boundary-layer thickness, has been exhaustively studied 
for some time, while the other regimes have not. These latter domains 
differ from the continuum flow both in the form of the equations of 
motion and in the boundary conditions. The change in the boundary con- 
ditions appears in the form of temperature and velocity discontinuities 
between a solid boundary and the fluid immediately adjacent to it. 

The slip regime may be loosely defined to include flows such that 
the ratio of mean free path to boundary-layer thickness (which is shown 
in reference 1 to be proportional to the Mach number divided by the 
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square root of the Reynolds number) is between 0.01 and 1. It may be 
noted that this parameter is the Knudsen number based on boundary-layer 
thickness. In this regime the flow may be defined by the Burnett equa- 
tions (references 2 and 3), which represent the second approximation to 
the Boltzman equation, the first approximation being the familiar 
Navier-Stokes and energy equations. The boundary conditions were first 
given to the corresponding degree of approximation in reference 4, where 
it is shown that, in spite of the fact that the Burnett equations are of 
higher order than the Navier-Stokes relations, the same number of bound- 
ary conditions are required. 

Several authors have discussed special problems in slip flow (ref- 
erences 4 to 6, for example) such as the Couette flow between flat plates 
and concentric cylinders and an incompressible boundary- layer and 
stagnation-point flow. In references 7 to 10, slip flow over a flat 
plate is discussed from another point of view by analogy with the 
Rayleigh problem of an infinite plate suddenly set into steady motion 
in a viscous fluid. This approach is relatively simple, but is limited 
by questions of the validity of the analogy for slip flow. 

In the present paper, the first-order solution of the compressible 
flat- plate boundary-layer problem is found for variable Prandtl number 
and arbitrary temperature-viscosity law. It is assumed only that the 
specific heats are constant. The solution is a perturbation on the 
known continuum solution. Although the Burnett equations are the appli- 
cable ones in the slip regime, these expressions are Just the f ami liar 
continuum, equations of motion with higher-order shear and heat- flux 

terms added. These added terms are at least of order 6 2 f e = ), 

, . V . VSe/ 

(All symbols are defined in appendix A.) The corresponding boundary . 
conditions (reference 4) contain terms of first order in 6. Hence, to 
find the first-order effect of slip, the continuum equations of motion 
may be used. These equations are subjected to the boundary- layer 
assumptions and expanded in powers of s to give two sets of equations, 
one for the zero- order and the other for the first-order quantities. 

The same procedure is followed with the boundary conditions. The zero- 
order system is then seen to be the usual continuum boundary- layer sys- 
tem. The first-order system is then solved explicitly in termB of the 
zero- order quantities. For the convenience of the reader who is prima- 
rily interested in the results of this investigation, the rather lengthy 
mathematical details of the analysis have been relegated to the appendi- 
ces. The investigation was conducted at the NACA Lewis laboratory. 
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AMLYSIS 

•* Differential Equations 

If it is assumed that there exists a thin boundary layer close to 
the wall in the fluid, and appropriate nondimensional coordinates are 
introduced (appendix B) , the equations of steady motion and the equation 
of state (herein assumed to he the general gas lav) reduce to the famil- 
iar boundary- layer equations, which are correct to order e (equa- 
tions (Bl) of appendix B). It may next be assumed that the velocities 
and thermodynamic properties of the fluid can be expanded in powers of 
the s m all, but nonzero, quantity s. Thus, for example, 

u = °u + s^u + . . . 

T = °T + s^T + . . . 

In order to find the differential equations governing these zero- 
and first-order quantities, equations (l) are substituted into the 
boundary- layer equations (Bl) and the coefficients of the zero and first 
powers of e are equated to zero. The resulting two sets of relations 
are equations (B2) and (B3) of appendix B. The first set contains only 
zero-order terms and is Identical with the usual boundary- layer rela- 
tions (equations (Bl)), while the other set governs the first-order 
terms; it is the pr imar y object of the present report to solve them Bub- 
* ject to the appropriate boundary conditions. 



Boundary Conditions 

The boundary conditions at the plate have been discussed in detail 
by Schamberg in reference 4 (see also reference 6) and are given as 
equations (B4) of appendix B. If these are treated in the same manner 
as the differential equations, that is. If they are made nondimensional 
and expanded in powers of £, equating the coefficients of the zero n.nrl 
first powers of e to zero, two sets of boundary conditions (equa- 
tions (B5) and (B6)) result. Finally, conditions far from the plate 
must be specified. At first glance, It might be expected that the con- 
ditions are simply that the properties approach the undisturbed stream 
values. To zero order, namely in the conventional boundary- layer 
theory, this is nearly true; but when higher-order approximations are 
desired, the mutual interaction of the boundary layer and outside stream 
must be considered. If, as is the case in this report, it is desired to 
use the known continuum boundary- layer solutions as the zero-order solu- 
tion to the present problem, it is found that the zero-order solution is 
v such that, as the outside of the boundary layer is approached, al 1 the 

fluid properties approach the undisturbed stream values except the ver- 

^ tical velocity. This component, which is of magnitude l/\J Re and is 
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hence negligible to zero order, is a manifestation of the fact that the 
boundary layer causes a curved surface to be presented as the effective 
wall of the relatively inviscid outside stream flow. This outer flow, 9 

which to the present appro xima tion satisfies the Prandtl-Glauert poten- 
tial equation, can be solved subject to the known boundary condition on 
the vertical velocity at the edge of the boundary layer (see appendix B, 
equations (B8) and (B9)). The resulting axial velocity, and hence the 
pressure and temperature (the outside stream remains isoenergetic to 
first order) at the edge of the boundary layer, may then be found (equa- 
tions (BIO) and (Bll)). These terms are of order e for supersonic 
flow and of order e*, and hence zero to the present approximation , for 
subsonic flow. It then follows that the appropriate first-order bound- 
ary condition is singly that the first-order fluid properties in the 
boundary layer must reach these values (equations (B13) and (B14)) at 
the edge of the boundary layer. Actually, the first-order solution will 
be seen to behave in such a man ner that these conditions yield insuffi- 
cient information to make the solution unique. Hence a suitable added 
condition must be introduced. A sufficient one is that the momentum and 
displacement thicknesses be finite. It will be Been that this condition 
effectively requires that the various fluid properties approach their 
stream values exponentially. 


Solution of Differential Equations 

The zero-order system • (equations (B2), (B5), and (B7 ) ) is the usual 
boundary- layer system and has been treated exhaustively (references 11 
to 13) under various assumptions regarding the Prandtl number and the 
viscosity-temperature relation, and its solution is partially described 
in appendix B. The first-order system (equations (B3), (B6), and (B13)) 
can be Bolved explicitly in terms of the known zero-order solution if 
similarity is again assumed in the usual form. The details of the 
rather involved analysis are given in appendix B following equa- 
tions (B18). The resulting solution is, with use of equations (B19) 
and (B3l) , 





MCA TN 2818 


5 


^(x.y) = % (°n(T v )ai. " f 4t ’ J' V F(t)dt ) 

***> -- js \ (°»Wi - f £ 4t - £ ? y(t)dt ) 

Mx ' ri ' ‘ ° £ (^* ^mS 1 )] 


°£ Pr dt 
2°\°I 


KoOCci - a x ) ^°T „(0) 


4 


°K Pr dr 
2°X°T 


(t - b) ( r - i)M 2 Pr°u+ ^3 - 

L ^nJ 


Also, it should he noted that 


K = 0 for subonsic flow 


and for supersonic flow 


K = - 


krr O 


M\/r(M^-i) 


iVr(M 2 -i) 


r — 8* 

V®^3f 


where 6* is the zero- order displacement thickness. 


From these expressions, the local skin friction and local heat 
transfer can be found. The result is, to first order (appendix C), 



6 


NACA TN 2818 


o f jBT x = °(c fA /R^) + 6 x l(c f ^R^} = [2 0 ii°u T) + e x (2M/^/r(M 2 -l) K 2 ]^ 



where , again, equation (3) holds. 


DISCUSSION OF RESULTS 


' In order to see where the various terms in equations (2) and (4) 
arise, it may be ob served that, if there were no slip, a^ and cq_ would 
be equal to zero (see the wall boundary conditions, equations (B4), (B5 ) , 
and (B6 ) ) , while if there were no interaction between the boundary layer 
and the stream outside it, K would be zero. Thus, for example, for 
subsonic flow with no slip the entire first-order solution vanishes. 


With reference to equations (4), it may be observed that, for sub- 
sonic flow, the second approximation (first-order solution) contributes 
nothing to the skin friction, while for supersonic flow this approxima- 
tion has a contribution. This skin- friction effect is not due to the 
slip at the wall at all but is rather simply due to the interaction 
between the boundary layer and the free stream. It therefore represents 
the effect of a slightly decelerating flow, generated in this case by 
the bo undar y layer itself, on the boundary layer. On the other hand, 
the change In heat transfer due to the first-order solution is entirely 
a slip effect and there is no basic difference between the subsonic and 
supersonic cases. 


Some comparison of these results with those of other authors may be 
made. For incompressible flow, the present solution for -hi(x,y) is 
the same as that found in reference 5, and both predict zero first-order 
shear (for subsonic flow). The Rayleigh problem (references 7 to 10) 
also yields solutions predicting zero first-order shear. The heat- 
transfer solution found In reference 10 is incorrect because the fric- 
tional dissipation which was neglected is the source of the thermal 
energy involved in the first-order heat transfer. The solution given in 
reference 8 for heat transfer is alBO incorrect because, in the notation 
of that report, as Te (adiabatic equilibrium temperature) is a func- 
tion of 9 (time), the solution for Tg, from which the heat tranfer 

was found, is in error. If, however, the analysis of reference 8 were 
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carried through, without error, the same result for first-order heat 
transfer as that found in this report would he obtained (appendix D). 

If it were desired to find only the second approximation to the 
laminar compressible boundary layer for constant wall temperature, 
assuming no slip (see, for comparison, reference 14, where this problem 
is considered for the case of incompressible flow), the solutions 
already obtained would apply except that the factors a^_ and c^_ 

appearing in the solutions (equations (2) and (4)) introduced by the 
slip conditions would in this case be set equal to zero. In this con- 
nection, it may be observed that, when, at the start of the analysis, 
the Navier-Stokes equations were expanded in powers of e, only even 
powers of e appeared. Also the wall boundary conditions involve no 
odd powers of e, and in fact are, except for the zero-order temperature, 
homogeneous. A superficial observation of the conditions to be applied 
at the outside of the boundary layer would perhaps lead one to expect 
that these also were, when expanded, independent of odd powers of 6 
(and would in fact be homogeneous except for the zero-order terms) and 
that therefore it would be proper to solve the problem by expanding ip. 
powers of e?. That this is not the case, at least for supersonic flow, 
has been shown here and the reason lies entirely in the conditions which 
arise in matching the boundary layer with the outside, essentially 
inviscid, stream. 


NUMERICAL RESULTS 

The first-order skin friction and heat transfer may be readily cal- 
culated provided the zero-order skin friction, heat transfer, and dis- 
placement thickness are known. The data of reference 13 have been used 
as the known zero-order solution because they are the best available 
continuum, solution of the flat-plate laminar boundary- layer problem for 
high-speed flow. The notation of reference 13 is related to the present 
te rmin ology by the following relations: 
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Present report 

Reference 13 

-2KM/^r(M 2 -l) 


T»* 


V 


V 

T i 

°(c£AyRe x ) 

(c f is/Si ) 8 

M 

M 

°( Nu x \ 

I Nu x \ 

\3s /i?r" */Re x / 

V$/Pr >\fRe)§ 


The results of the calculation are presented in table I, which is 
to he used with equations (4), together with, for comparison, the zero- 
order solutions as given in reference 13. For computational purposes, 
a-j^ and Cj_ were taken as 1.253 and 2.870, respectively, corresponding 

to values of the Maxwell reflection coefficient of 1.0 and of the accom- 
modation coefficient of 0.9 (reference 6). Computations are shown for 
an ambient stream temperature T^* of 400° R, which is fairly typical 

of high- altitude flight. The quantity e x is, in free air, for a 
standard atmosphere (reference 15), given by 



where H is the altitude. The variation of a(H) with altitude is 
shown in figure 1. Thus, for example, at an altitude of 200,000 feet, 
if the Mach number is 10, the plate temperature is 1200° R, and 
x* = 1 foot, the first-order solution shows a 10-percent increase in 
the shear and the same amount of decrease in the heat transfer. At 
300,000 feet and the same Mach number, temperature, and length, the 
first-order solution indicates 60-percent increase in shear and a simi- 
lar decrease in heat transfer. Of course, in this latter case, e x is 

nearly one half and is certainly too large for the theory presented 
herein to give more than a qualitative answer. 
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CONCLUDING REMARKS 

The first-order boundary- layer solution for a flat plate at con- 
stant t emp erature has been calculated for high-speed motion in a semi- 
rarefied gas. There are actually two effects involved. One is the 
effect of slip itself, as reflected in the boundary conditions at the 
wall, and the other is a result of the interaction between the boundary 
layer and the stream outside it. These two effects are independent of 
each other and the effects in the latter case are not the same in the 
subsonic and supersonic regimes. The slip causes a decrease in the heat 
transfer at the wall but has no effect on the shear to first order, while 
the second effect causes, for supersonic flow, an increase in the shear 
and does not affect the heat transfer. For subsonic flow, the second 
effect causes no change in either the shear or heat transfer to first 
order. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 6, 1952 
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APPENDIX A 
SYMBOLS 


The following symbols sire used in this report: 
AjA^A^Aj constants of integration 


a l 


b 


c f 


C 1 


F,f 

G 

H 

K 

L 

l 

l * 


M 


nr 2-a , 

V 2 IT’ where 

(reference 6) 
lim °v(x,T)) 


a 


T) -*■•» 


is Maxwell's reflection coefficient 


constants 

local skin-friction coefficient 

« 

specific heat at constant pressure 
1 

where c is the accommodation coef- 
ficient (reference 6) 

functions defined by equations (B26) and (B3l) 
function defined by equation (D4) 



altitude, ft 


constant, * fo — lim *1* °v(x,y) 
\|M -1 y-Ko 


characteristic length, say, length of plate 


l*/ L 

quantity defined by equation (D5)j approximately the molec- 
ular mean free path 

Mach number of undisturbed stream 


m 


m*/L 


r 


4 


Ir 
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m* - 

Nu, 

P 

P* 

Pr 

R 

Re 


Re x 

r,s,t 

S 

T 

T a 

T v 

rp* 




T* 

w 


u,v 

u*,v* 


quantity defined by equation (D5) 
local Nus selt number 
nondimensional pressure, P*/p£ 
pressure 

Prandtl number in undisturbed stream 
gas constant 

Reynolds number based on characteristic length and undis- 

p*u*,L 

turbed stream values, — 


local Reynolds number along plate, 


p£u*x* 


dummy variables 

function defined by equation (D4) 
nondimensional temperature, T*/T* 

*j/t; 

temperature 

continuum adiabatic wall temperature 
specified wall temperature 


nondimensional cartesian velocities; u = v = 


u^ 


v* 

u* 


*yse 


cartesian velocities, u* being parallel to the plate and 
v* normal to the plate 

nondimensional cartesian coordinates; x = x*/L, y = -^ ry/Re 

cartesian coordinates, x* being parallel to the plate 
y* normal to the plate 


x*,y* 
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function of altitude defined "by equation (5) and 
reference 15, ftV^ 

function defined by equations (B19) and (B3l) 
ratio of specific heats 
displacement thickness 

mVt 

slip parameter, 
time 

y/VT 

momentum thickness 

nondimensional heat conductivity, 1*/^* 



heat conductivity 
nondimensional viscosity, M-*/k£ 

\ dT } rj—Ogi 
\ dT^ / ij_Oiji 


viscosity 

function defined by equations (B16), (B19), and (B3l) 

nondimensional. density, p*/ P* 

density 

function defined by equations (B19) and (B3l) 


shear 
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§,<?> velocity potentials 

Subscripts: 

o» undisturbed stream 

s stream conditions outside boundary layer 

x,y,f;*,T} partial derivatives with respect to respective coordinate 

except for Re x , e x , and ITu^., which are defined elsewhere 

Superscript: 

* physical quantity 

Pre superscript : 


0,1 


zero and first approximations, respectively 
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APPENDIX B 


ANALYSIS 

The Differential Equations and Boundary Conditions 

On application of the houndary-layer assumptions and introduction 
of nondimens ional coordinates, the equations of steady motion and the 
general gas law reduce to these familiar houndary-layer equations , cor- 
rect to order e , 


(pu) x + (pv) y = 0 


P x + pCUU^+VUy) - (jlUyJy = 0 


p 


y 


= 0 


> (Bl) 


r-i 


uP x - p(uT +vT ) + — (XT ) 
x y' Pr y 


(r-l)M 2 pu y 2 = 0 


P = pT 


J 


In the transformations involved the coordinate and velocity normal to 
the wall are multiplied hy l/ -/Re . Thus the physical and velocity 
fields are stretched in one direction, namely, across the houndary layer. 
ThiB Is, indeed, ,)ust the usual houndary- layer procedure. If equa- 
tions (l) are substituted into equations (Bl), the resulting expres- 
sions are expanded, and the coefficients of the zero and first powers 
of e equated to zero, two sets of equations result. These are 


/O O \ /O O \ _ 

( P u) x + ( p v) y = 0 


1 °P X + 0 b (°uV° t V - <V>vV - 0 


rM 2 


P = 0 

y 


VB2) 


°u°P x - °p( 0 u 0 T x + 0 v°T y ) + i (°X°T y ) y + (r-l)M 2 °|i°u y 2 = 0 


°P = °p°T 


¥ 
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(°p 1 a+ 1 p°u) x + (°p 1 v+ 1 p 0 v)y =0 " 

1 1_ , ° /° 1 1 o o 1 ,1 o ^ ,1 ,o o o o y ,o 1 ,1 o v _ 

2 T* + p( U U x + U OJJ.+ V Uy+ V Uy) + p( U U x + V Uy) - ( P Uy+ \1 Uy)y = 0 

rtr 


p y= 0 


~ ( 1 u°P x +°u 1 P x ) - °p ( 1 u°T x + 0 u 1 T x + 1 v 0 Ty+°v 1 Ty) - X p (°u°T x +Vt v ) + 




x * -y' 


P? ( 0xLT y +1 ^° T y)y + (r-l}M 2 ( 1 p°Uy 2 +2 0 p°Uy 1 Uy) = 0 
h? = °p 1 T + 1 p°T 


(B3) 


The boundary conditions at the plate are (references 4 and 6 ) , to 
first order. 


v*(x*,0) = 0 




(B4) 


T*(x* 




where the terms in brackets are evaluated in terms of conditions in the 
fluid adjacent to the wall. 


If these are written in terms of the dimensionless quantities, there 
follows, correct to order e, 


v(x,0) = 0 

u(x,0) =60 ax-v/OV u y ) 

T(x,0) = T w +60 c-l-/^ T^ 

These relations may next be expanded in powers of e to give the two 
sets of boundary conditions at the wall. These are 
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°v(x,0) = 0 
°u(x,0) = 0 
°T(x,0) = T w 

^(x^O) = 0 

■^(xjO) = 3^ a L °Uy 

, o 

1,1 U o 

- SJ ®1 T y 

Equations (B5) are the zero -order or continuum boundary conditions, 
while equations (B6) govern the first-order quantities. 

Finally, to complete the definition of the problem, conditions far 
from the plate must be specified. Because it is desired that the zero- 
order problem be the usual continuum one, and because this introduces 
no inconsistency, the zero -order boundary conditions may be immediately 
specified. These are 

°u(x,"} - °P(x,<*>) = °T(x, w ) = 1 (B7) 

If it is now noted that the zero -order system (completely specified 
by equations (B2), (B5), and (B7 ) ) is actually the desired familiar con- 
tinuum one, it is then known that 

°v(x, °°) = £ -/x (B8) 

where b is a known constant. The outside stream satisfies the Prandtl- 
Glauert relation 



(l-M 2 ) Vx* + Vy* = ° (B9) 

where <j> is the perturbation (from the undisturbed stream) velocity 
potential. The solution of equation (B9) subject to condition (B8) 
applied at the edge of the boundary layer will be found first for super- 
sonic flow. The applicab le gen eral solution of equation (B9) is that <P 
is a function of (x# - y*) • Thus 
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Then, with use of equations 
■boundary layer 


1 


u 


s 



(1), (B7), and (B8), 

_ b 

2MVr(M 2 -l)x 


0(e)] 

at the edge of the 


(BIO) 


Consideration will now be given to subsonic flow. A solution of 
equation (B9) corresponding to a linear source distribution of strength 

varying as l/ yields 


u* = - b 


,yg a /A* 2 - 0 *-^ ' * 
\l x*Z- ( M 2-l)y*2 


+ i£ 


_lL_ . * - 

V x* 2 - ( M 2-l)y*2 


or, in nondimens ional coordinates, 
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6v. 


X 


M^rd-M^y 




This last relation is of the same form as equation (B8) and hence 
these two expressions may be considered to define the desired flow out- 
side the boundary layer for subsonic flow. Thus, by the preceding 
equations, b-j^ is of order e. Although b]_ Is proportional to e/M, 
which may apparently be large if M approaches zero, It is seen from 
the definition of e that b^ is proportional to l/-/Re, which must 
be small in order that the boundary-layer assumption will hold in the 
first place. Hence, Ug - 1 is of order e 2 . Therefore, for subsonic 
flow, 

^ = 0 (B12) 

The pressure and temperature perturbations will now be described. 
For Isentropic flow, which condition is satisfied to the present approxi- 
mation by the outside stream, 


T s - 1 = - (r-l)M 2 (u s -l) 


P s - 1 = - TM 2 (u s -1) 


Thus, for the first-order boundary conditions at the edge of the 
boundary layer, 

1 , S K 

X u(x,-) = —= 

Vx 


W) - - 

Vx 


1 p(x,<») = - 


Vx 


(B13) 


where 
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K = 0 for subsonic flow 


K = 


-1 


2MVr(M 2 -l) M-/r(M 2 -l) 


[Vx °v(x,»}] 


(B14) 

for supersonic 
flow 


In practical calculations, when the zero-order solution is known, 
K may be found in the following way: 


K = - 


1 

M-/r(M 2 -i) 



V s ; as* 


where 6* is the continuum displacement thickness and is known to "vary 
as Vx* Hence, 


K a - 




2MVr(M 2 -l) 


5 


(B15) 


Solution of Differential Equations 


First, the zero-order system (equations (B2), (B5) , and (B7)) may 
be considered. From the third of equations (B2) and the boundary con- 
dition equation (B7), °P(x,y) = 1. How, similarity may be assumed in 
the usual manne r and there can be written, in the manner of Emmons and 
Bralnerd (reference ll) : 


\ 

°u(x,y) = °u(t]) 

°v(x,y) = |ji°u(ti) - °£(r))j \ 

°T(x,y) = 0 T(tj) 

J 


(B16) 


where 


t) = y/Vx. 



20 


NACA TN 2818 


Then the system becomes the familiar one 


°T \°T/ 


°*\ + 2°T(°A n ) n = 0 
1 + h + S^-ljM 2 Vu/ = 0 


Vt, 


°T Pr 


VTJ 


J 


(B17) 


°u(0) = °«(0) = 0 

°T(0) = T^ 

°u(-) = °T(-} = 1 


J 


(B18) 


This system or its equivalent has been solved by a number of 
investigators under various assumptions regarding the Prandtl number 
and the temperature-viscosity lav (see, for example, references 11 to 
13) and, hence, need not be considered further. 

Next, the first-order system (equations (B3), (B6), and (EL3)) may 
be considered. From the third of equations (B3) and the third of equa- 
tions (B13) it follovs Immediately that 


W) - 

Vx 


Next, analogously vith equations (B16), and following Karman’s 
suggestion (reference 16), it may be assumed that the first-order 
quantities have the form 


1 u(x,y) = Mn) 

Vx 

^(x,y) = ^ [*(,) 
1 T(x,y) = -7= 1 P(t 1 ) 

■yx 



■\ 

J 


(B19) 
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and, hence , from the last of equations (B3), 


^(x.y) i 

T(n) 

Also, since the viscosity p and heat conductivity X are func- 
tions of temperature alone, ^p and ^X may he found by expanding p 
and X in a Taylor series. Thus, 


KrM 2 + 

°T(n) 


p(°T+e?-T) = p(°t) 



+ e^p 


Therefore, 


Similarly, 


where 


V(x,y) = vy 


1 X(x,y) = 




J 


(B20) 



Now the expressions given by equations (B19) and (B20) may be put into 
the appropriate equations (B3), (B6), and (B13) to yield this system 
for the first-order terms. 
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K = 0 


"5 f + + ^ (ktm 2 + + 2(v°u, + °l*\) t 

^ + "V 1 + " ° 5 + S] + & + (r4) w 

2(r-l)M 2 ( 1 P|J T 0 u 1) 2 + 2 0 ii 0 u I1 1 <l I1 ) = 0 


> 


+ 


1 C(0) = 0 

Mo) = V-lVt^ M^ 0 ) 
Mo) = ViV^ \(0) 
M«) = K 
M») = - (r-il^K 


(B21) 




> (B22) 


J 


Equations (B2l) may be Integrated explicitly in terms of the zero- 
order solution. From the first of equations (B2l) and the boundary con- 
ditions on °£(t)) and ^£(t]) (equations (B18) and (B22 ) ) , there follows 
immediately 


Mt)) 





(B23) 


With this result the last two equations of equations (B2l) can be 
integrated once to give 


O - •. 

cr + 2( 1 p t i T °u T) + V^) = 2A + Kq 

^ 1 P + ^ (^p)^ + 2(r-l)M 2 Vu^ff = ^ - 2(r-l)M 2 A°u - (r-l)M 2 K°uriJ 

(B24) 
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It can. readily "be verified that a solution of the homogeneous equations 
corresponding to equations (B24) is 

1 0(Tl) = °u^(ti) 


Vq) = VCtj) 


and, hence, in general, 




r n 


1 cr(n) - % 

Ai + 

f (t)dt 


1 


0 

\ 





1 . . o 

p(n) = 

Ai + 

f (t)dt 

+ F(q) 


L vJ 

0 

/ 


(B25) 


If these expressions are put into equations (B24=) , the result is 

(V + t^) 0 ^ - A + ^ 


I; + S 1 ) F + <°V +F rP - ?f - (r-i)M 2 °u(a + n) 


These equations may be solved to give 


^ 


I 


SPr at 
2°X°T 




+ f *2 - ^ * f) (r-D^n . 


pt 

°£Pr dr 
10 2°X°T 


At \ 


(B 26 ) 


and 


f(q) = 


A+Kn/2 


0 0 o 

* u q ^ 
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It may be noted that, in the commonly assumed case of constant 
Prandtl number, °p = °X and hence 

pn 

^Pr dt 



2 X T 


(Vy Fr 

(Vu n )^o 


In order to apply the boundary conditions for large tj, it is neces- 
sary to know the behavior of and P(t}) for large T). It is 

shown in appendix E (the result is derived for only cr ( rj ) while a simi- 

lar analysis holds for 1 p(r))} that for large 


1 o’(ti) = K 



- — — + exponentially 

decaying terms 





exponentially decaying terms 


(B27) 


and in the l i mit 


M-) = K 

1 P(«) = - (r-l)M 2 K 

Thus for all values of the constants of integration, the last two bound- 
ary conditions of equations (B22) are satisfied. Now the displacement 
and momentum thicknesses, 6* and 0, respectively, may be considered. 
These are, in nondimensional coordinates where the subscript s refers 
to stream values , 
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&*= (p s u -pu)dy = Vx I (p s u -pu)dT) 


6 = j pu(u s -u)dy = -yfx. / pu(u s -u)df] 


JO 


The contribution of the first-order slip term is 


» 

= V* J L P s °u s + °P s 1 u s - ( 1 p 0 u+°p 1 u)j 


AS' 




dt) 


A0 = s/x 


P 

J jj 1 p 0 u+°p 1 u) (°Ug — °u) + ^^(^g-^-ujJdT] 


> (B28) 


In evaluating these integrals , equations (B27) can he used because all 
that is desired is to find the conditions for them to be finite. The 
exponentially vanishing terms can contribute only finite amounts to the 
integrals. It then follows (see appendix F) that In order for AS* and 
AS to remain finite, 


h = 0 



(B29) 


Finally, the conditions at the wall (equations (B22)) must be applied. 
From these there follows 


A 1 = V T v) a i 

A 3 = VT w )( Cl - ai )V^ %(0) 


(B30) 


If equations (B26), (B29), and (B30) are put into equation (B25), there 
is obtained, finally. 
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These four equations constitute the solution of the first-order boundary- 
layer problem for slip flow when the wall has a constant specified tem- 
perature 
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APPENDIX C 


SKIN FRICTION AND BEAT TRANSFER 

The local skin friction and local heat transfer at the wall are 
given by 

V Be x c f = V Re x ^ 2 = 2 ^ u H^ti=0 

- 2 [ys + + 

W a/ 5 ^ - [ x *($-iS) (sp)L=o (^5=) ■ 5^$; 

■ [° x % + + ^X)]^ 

These become, with use of equations (B24) and the boundary conditions, 
op ^/Ri^= °(c f < V /R^) + e x 1 (e f - [f^\ + - 1) K 2 ), n=0 




v 5 ^" vv^y nv^y T *"M 


° X % + 6 X 


^ ai (r - lJl^PrCV ^) 2 
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APPENDIX D 


SOLUTION OF TEE RAYLEIGH PROBLEM FOR FIRST-ORDER SLIP 

The mathematical statement of the Rayleigh problem for slip in an 
incompressible fluid is (reference 8) (in the present statement of the 
problem it is assumed, as in reference 9, that the fluid rather than the 
plate is suddenly accelerated) 




u£* 


u 


y¥ 


= 0 


A 


u*(£*,0) - l*(uf*) y=0 

U*(£*>) = 

u*(0,y*) = u* 




J 


t£* 


m* P? 

£ y*y* ~ p* c 


P*Cp 


(Uy *) 2 




T*(£*,o) = Tj + **( T £*) y=Q y 

**(£*,-) = T* 

T*(0,y*} = T* j 


(Dl) 


(D2) 


Equations (Dl) can be solved exactly. If is chosen Buch that 

= x*/0.346u*, in order that the zero-order shear match the Blasius 

solution, there results, in nondimensional coordinates, 


u(x,y) = erf 


[ 0.294y | 


U/Re 0.346l 2 Re 
+ e erfc 


- 

0. 


294y 


Vx 


0.588i^/Re 


e 


from which there can be obtained 

c f <^Re^ = 0.664 + 0(«£) 


which is the result obtained in references 8 and 9, and which agrees 
with equation (4) of this report for subsonic flow. 
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Next, equations (D2) can be solved exactly "by standard methods. 
The solution is, in nondimensional coordinates. 


T(x,y) = T^ - S(x 


>y) + J* s(t, 


0)G(x-t,y) dt 


(D3) 


where 


*\ 


sfx.v'l = rp _ n + .(r-l)M 2 Pr2 
>} V 2^0.346) 


A* 

ds 

UO \J 




0.546Pr(r-y) 2 
2 4(x-s) 

e dr 


a/x-b 


and 


. _ 0 . 546Pry 2 


G<*,y) = — — — „ 


4x 


(; 


2 888 0.346m 2 RePr ; 


m*s/RePrx 

) 


m 2 RePr 


erfc 


£ 


Vx 


294y ft/ir + 

0 /x 0.588m^RePr, 


:) 


(D4) 




From this, to order Z (Z^/Re is of the same order as a), 
^|= " T^; L ( 1 _T v +B ) (°- 332 ) ^ - (0.332) 2 (r-l)M 2 ^ vs] 

However, if equations (Dl) and (B4) are compared, there follows: 


Z* , p* a-i^RT* 
L " 1 ~ p* L 


Re 


m* u* 

— = III = -Hr 


* c-j^RT* ec lV^w 


(D5) 


V5i 


Then 
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where 

£ , cos-1 Jf 

Knfz/^X-1 

The zero-order portion of equation (D6) has "been essentially found 
hy Emmons (reference 17), who has pointed out that, except for Pr » 1, 
this solution is not quite the same as that found from the usual 
boundary- layer theory. The first-order term in equation (D6) , however, 
exactly matches the corresponding term of equation (4) of the present 
report. 
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APPENDIX E 

BEHAVIOR OE ^(t}) FOR LARGE 7] 
From equations (B25) and (B26), 


Mn) - \ 






^Pr dr 
2°X°T 


e 


A 3 



First, It is well-known that the behavior of the zero-order solution 
is such that, as T) -*■ 


°u(t)) +1 

-biT] 2 

Ou(t]) - l -*-0 as e (b2_>0) 


o ' blTl2 

u u T j _> 0 as e 
°T(rj) ->1 

%(il) -»-0 as (%}) 


(E2) 


- b 

Most of the terms in equation (El) will be shown to vanish exponen- 
tially. For the terms multiplied by A^ and Aj, this follows immedi- 
ately by the use of equation (E2). Some of the other terms will be seen 
to vanish exponentially if equation (El) is multiplied by e and the 
product is shown to v ani sh. The terms in the last integral may be con- 
sidered. Altho ugh e^ °u^('q) -*0, the integral is infinite and hence the 

product is an indeterminate form, and it is appropriate to apply 
1* Hospital's rule. Thus, 
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Again, the integral approaches infinity and its coefficient is zero. 
Thus 



By the use of equations (E2), it is seen that the denominator 

approaches ~ - 1^ plus smaller terms. The numerator contains 

only terms of the form e^ °u r j, e^ iqe 1 ! °u^, and rje 7 ! all of 

which vanish. Hence, finally, 1-^ = 0 and therefore the last term in 

equation (El) vanishes exponentially. Next, the term involving K in 
the first integral is considered. This can be written 
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The last integral can he shewn to vanish exponentially. 

n 7 ) 

(VS) Bl 

Ui 


I 3 = lim 


(^S) n A . 

(°n°u,) 2 V ) 


0 


dt 


= lim 

T ] - * 00 


0 ^ 0 U T 1 ^ _ 

°z/2°\i°t - ±y °u°u T1 A 


°T(n-h) 


But, by equation (E2) , 


Hm I 3£t) 

^\%/20v.°T - 1 


Therefore, 


I 3 = lim 


[f (I - (T| - b) ) 


This is again an indeterminate form. Therefore, 


I 3 = lim 


Then, when these results are used and the indicated integration in equa- 
tion (E3) is carried out, it follows that, for large values of T), 


e]l_ 

fa A 

1-T) 

l°T 7 


'’'ffOl} = K +^A + vr j °Ufj 


P 7 ! 


dt 


°ti°u^ 


+ exponentially vanishing terms 


Finally, it can be shown that as tj approaches infinity, the term 
in this expression involving the integral vanishes as l/t). Hence, 

1 a(") = K 
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APPENDIX F 


BOUNDARY CONDITIONS AT INFINITY 

Consideration will now "be given to the expression for A 6 in equa- 
tion (B28) . From the form of the solution, the term (■ L p°u+°p'Ki) is 
bounded in absolute value hy some constant, bg • Then, with the knowledge 
that °u a * 1 and that °u - °u a is always nonnegative. 


(°p 1 u+ 1 p°u) (°u a -°u)dT} 


<b 2 (l-°u)dr) 

Jo 


which is finite. The remaining term In A 9 is, if equations (B19) and 
(B27) are used, 


1 4 = V* 


(**¥) f 4°\f 


dt 


(FI) 


But 


lim vq 


o r \ 

u o ds | 

o u n I oo r 

T do [i u„ / 


and is always such that, for some bj, 


u o 
11 5T “1J 


ds 


>b, >0 if ti>1 

T ' ^0 p u_ 


Hence , 


A ® 

* v4 + ¥) b 4 £ =v4 + “) 


b, lim (in tj) 


Therefore, A 9 is only bounded if 


A = 


bK 

T 


(F2) 
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Now, if equation (F2) is used and equations (B27) are put into the expres- 
sion for AS* (equation (B28)), the result is 



The first integral is finite because its integrand approaches zero 
exponentially. If the same procedure is observed as in the case of the 
expression in equation (Fl), it follows that the second integral is 
infinite and, therefore, Ag must be zero in order that AS* be finite. 
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TABLE I - VARIATION OF SKIN FRICTION AND HEAT TRANSFER 
WITH WAIL TEMPERATURE AND MACH NUMBER FOR NOMINAL 


STREAM STATIC TEMPERATURE OF 400° R 


(°R) 

M 

°(c fA /Re^) 

1 ( c f V^®x) 

°(Nu/VRe^) 

1 (Nu/VRe^) 

1182 

4.55 

0.591 

0.958 

0.311 

0.881 


6.31 

.572 

.738 

.319 

.525 


9.70 

.538 

.684 

.325 

.417 


14.83 

.481 

.752 

.313 

.321 

2246 

6.68 

0.531 

1.460 

0.312 

1.356 


10.20 

.503 

1.085 

.319 

.640 


15.50 

.461 

0.967 

.305 

.450 




a(H) 
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Figure 1. - Variation with altitude of a(H) for use in 

equation (5) . 


NACA - Langley Field, Va. 


